Integrable Quantum Computation 
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Abstract Integrable quantum computation is defined as quantum computing 
via the integrable condition, in which two-qubit gates are either nontrivial 
unitary solutions of the Yang-Baxter equation or the Swap gate (permutation). 
To make the definition clear, in this article, we explore the physics underlying 
the quantum circuit model, and then present a unified description on both 
quantum computing via the Bethe ansatz and quantum computing via the 
Yang-Baxter equation. 
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1 Physics underlying the quantum circuit model 

Computers are physical objects, and computations are physical processes, which 
may be one of Deutsch's famous quotes [1]. This statement gives rise to a natural 
question: what is the physics underlying the quantum circuit model? Deutsch's 
answer may be that quantum computing votes for the many-worlds interpretation 
of quantum mechanics. Nielsen and Chuang's answer to it stays the same in the 
two editions of their popular textbook [2]: "outside the scope of the present 
discussion". In our current understanding [3], the quantum circuit model can 
be viewed as a generalization of the multi-qubit factorisation scattering model in 
integrable systems [3], and the scalable quantum computer is thus supposed to 
be an exactly solvable model satisfying the integrable condition [4]. 

It becomes well known in quantum information science after DiVincenzo's 
work [5] that an arbitrary A^-qubit quantum gate can be expressed as a sequence 
of products of some two-qubit gates. Hence a quantum circuit is described as 
a network of two-qubit gates, which may be called the locality principle of the 
quantum circuit, see Preskill's online lecture notes j6j. On the other hand, it has 
been widely accepted for a long time in integrable systems |4j that many-body 
factorisable scattering can be expressed as a sequence of two-body scattering. 
As two-qubit quantum gates are considered as two-qubit scattering operators, 
the quantum circuit is thought of as the generalized factorisable scattering model 
defined as a generalization of the factorisation scattering model [1]. 

Sutherland in his book [4] describes three types of the integrable conditions 
which are capable of yielding exactly solvable quantum many-body systems: the 
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first is the Bethe ansatz; the second is the Yang-Baxter equation; and the third 
is the quantum Lax equation. In physics, hence, we define integrable quantum 
computation as quantum computing via the integrable condition so that we can 
obtain a unified description on both quantum computing via the Bethe ansatz [3] 
and quantum computing via the Yang-Baxter equation [7j. Note that we do not 
deal with quantum computing via the quantum Lax equation in this paper. 

Section 3 is a prehminary introduction on quantum computing via the Bethe 
ansatz [3]. The main reference is Gu and Yang's paper [8j on the application 
of the Bethe ansatz to the Fermion problem with delta-functional potential. 
Given a quantum many-qubit system, i.e., the Hamiltonian is known, if the Bethe 
ansatz is satisfied, then it is an integrable model to yield many-qubit factorisable 
scattering as a sequence products of two-qubit scattering matrices in which two- 
qubit gates satisfy the Yang-Baxter equation. 

Section 4 is a brief sketch on quantum computing via the Yang-Baxter equa- 
tion [7]. Given a nontrivial unitary solution of the Yang-Baxter equation, the 
Hamiltonian of an integrable system can be constructed in principle. Addition- 
ally, the paper [7] is written in the way to emphasize how to derive unitary 
solutions of the Yang-Baxter equation from unitary representations of the braid 
group. In view of the research [3], however, it has been distinct that integrable 
quantum computation is a research subject independent of quantum computing 
via unitary braid representations. Hence this section is a refinement of the re- 
search [7] from the viewpoint of integrable quantum computation. 

In mathematics, integrable quantum computation specifies a type of quantum 
circuit model of computation in which two-qubit gates are either the Swap gate 
(permutation) or nontrivial unitary solutions of the Yang-Baxter equation. 

2 Quantum computing via the Bethe ansatz 

Let us consider the model of qubits (spin-1/2 particles) in one dimension 
interacting via the delta-function potential, with the Hamiltonian given by 



where c > (c < 0) means repulsive (attraction) interaction, see [8j for the 
detail. Interested readers are also invited to refer to Bose and Korepin's article 
[9] for the physical reason why this model is interesting in experimental quantum 
information and computation. 

For two qubits at xi and X2, the outgoing wave in the region xi < X2 with 
respective momenta ki < k2 is given by 
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the incoming wave in the region xi < X2 is given by 




amc = "21, 



(3) 



2 



so that the two-qubit gate (i.e., the scattering matrix) Ri2{k2,ki) defined by 
ctout = Ri20iinc has the form 

n (I. h\ ^(^2 - fci)(-fi2) +c 

where the permutation operator P12 exchanges the spins of two flying qubits. 

For three qubits at xi, X2 and X3, the outgoing wave in the region xi < X2 < X3 
with respective momenta ki < k2 < k^ is given by 

aout e^('=i-i+'=2-2+^3X3)^ = «^23, (5) 

the incoming wave in the region xi < X2 < X3 is given by 

SO that Oout is determined by Omc in the way 

Oout = Rl2{k2, ki) ^23(^3, ^1) Ruik-s, k2) Oinc, (7) 

or in the other way 

Oout = -^23(^3, ^2) R23{k3, ki) R23{k2, ki) Oinc, (8) 

which give rises to the consistency condition (the Yang-Baxter equation) on the 
three-qubit scattering, 

Rl2{u) R23{U + V) Ru{v) = R23{v) Rl2{u + v) R23{u) (9) 

with u = ki — k2, V = k2 — k^ and u + v = ki — k^. 

The outgoing wave of two qubits at xi and X2 in the region xi < X2 with 
respective momenta ki < k2 is given by ([2]), but the incoming wave in the region 
X2 < xi is given by 

a,„,e*(^^^^+'=^^2), a,„e = (-A2)a2i, (10) 
so that the scattering matrix of two qubits in different regions has the form 

aout = Rl2{k2, ki){-Pi2)ainc- (11) 

Similarly, the outgoing wave of three qubits at xi, X2 and X3 in the region xi < 
X2 < X3 with respective momenta ki < k2 < k^ is given by (l5|), but the incoming 
wave in the region X3 < X2 < xi has the form 



"321 e 



,(k,.,+k,.,+H.,)^ a.„c = (-Pi3)«32i, (12) 



where the permutation operator P13 exchanges the spins of the first qubit and 
the third one, so that the scattering matrix of the three qubits in different regions 
is given by 

aout = Ri2{k2, ki) ^23(^3, ki) Ri2{kz, k2) {-P12) {-P23) (--Pi2)amc (13) 
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where P13 = Pu P23 P12 is used. 

In term of the notation R{u), the identity operator Id and the permutation 
operator P, the two-body qubit gates R12 and R23 have a new form 

R(u)= ~^ (c + iuP), Ri2 = R(E>Id, R23 = Id(S>R, (14) 
lU + c 

where the two-qubit gate R{u) satisfies the unitarity condition, 

R{u) R\u) = R\u) R{u) = Id, R\u) = R{-u). (15) 

Hence the multi-qubit factorisable scattering with the delta-function interaction 
has been recognized as a quantum circuit model consisting of the Swap gate P 
and the two-qubit gate R{u) (fH|) . 

Besides the form (I14p . the two-qubit gate R{u) has the other interesting for- 
malism 

i?(^) = -e-2^^P_-P+, P± = ^^- (16) 

with If = arctan ^ and —^<<p<0. AtM = — c, the two-qubit gate R{u) is the 
y/Swap gate given by 

R{-j) = -y/Swap, ^JSwap = P+ + iP- (17) 

which is an entangling two-qubit gate but can not yield all unitary 4x4 matrices 
by itself (even with the Swap gate) p^ . 

Furthermore, the two-qubit gate R{'^) has the form of the time evolutional 
operator U (if) of the Heisenberg interaction 5i • S2 between two qubits 

R(^) = _e-*f C7(-2(^), U{^) = (jg) 

modulo a global phase factor with the permutation P = 2Si • 5*2 + ^ so that R{(p) 
and U{ip) are two equivalent two-qubit gates in the quantum circuit model. 

As is realized by DiVincenzo et al. [11], universal quantum computation can 
be set up only using the Heisenberg interaction when a logical qubit is encoded 
as a two-dimensional subspace of eight-dimensional Hilbert space of three qubits. 
Hence, the quantum circuit model in terms of the Swap gate P and the two-qubit 
gate R{u) (I14p is also able to perform universal quantum computation, if and only 
if a logical qubit is chosen in a suitable way. 

3 Quantum Computing via the Yang-Baxter equation 

A quantum circuit model in terms of both a nontrivial unitary solution R{u) 
of the Yang-Baxter equation ([9]) and the Swap gate P is called the generalized 
factorisable scattering model in the paper. Once it is given, the problem becomes 
how to find out the underlying physical system, which is thought of as the inverse 
problem of quantum computing via the Bethe ansatz. 
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In the literature, the Yang-Baxter equation has the usual formalism, 



{Rix)(g)l2) (l2«>^(xy)) {R{y)(E)l2) = il2^R{y)) iR{xy)(g)l2) il2(S)R{x)) (19) 

where the two-qubit gate R{x) is a linear operator on the Hilbert space of two 
qubits, i.e., R :% ^ % (^T-L, IL2 denotes the 2x2 identity matrix, and x, y 
are called the spectral parameter. Take x = e" and y = e" , then the formalism of 
the Yang-Baxter equation ([9]) can be derived. As a two-qubit gate, the solution 
R{x) of the Yang-Baxter equation (|19p has to satisfy the unitary condition 

R{x) i?t (x) = ijt (2;) R[x) = pl4 (20) 

with the normalization factor p. 

The construction of the Hamiltonian via a nontrivial unitary solution of the 
Yang-Baxter equation (jl9p is a very flexible or subtle process with physical rea- 
soning, namely, there does not exist a universal law to guide such a construction, 
partly because the same solution of the Yang-Baxter equation can be yielded by 
many different kinds of physical interactions. 

The two-qubit gate , as a rational solution of the Yang-Baxter equation 
(I19p . has the form 

1 u 

R(a) = (IL4 + aP), a = i- (21) 

1 + a c 

which is a linear combination of the identity 1 and the permutation P. In terms 
of the variable ip, the two-qubit gate (I14p can be expressed as an exponential 
formalism of the permutation P given by 

R[ip) = e'^^-^^e'^^, tan(/7 = -, (22) 

c 

hence the Hamiltonian for the generalized factorisable scattering model can be 
chosen as the Heisenberg interaction Si ■ S2 between two qubits. On the other 
hand, it is easy to examine the delta-function interaction ([T]) as a suitable Hamil- 
tonian underlying the two-qubit gate (|14p . but it is not so explicit to realize the 
delta- function interaction ([1]) from the formalism of the two-qubit gate (jl4p . 

In [7], an approach for the construction of the Hamiltonian is presented via a 
two-qubit solution R{x) of the Yang-Baxter equation ()19p . With an initial state 
ip, the evolution state ip{x) determined by R{x) is given by ip{x) = p^'i R{x)ip. 
The Shrodinger equation has the form 

i^^ = Hix)i;{x) (23) 
with the Hamiltonian H{x) given by 

H(x) = i—{p-^R)p-^R-Hx) (24) 
ox 
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which is time-dependent of x. 

In terms of an appropriate parameter instead of x, the Hamiltonian H{x) ()27p 
may be transformed to a time-independent formahsm. For example, the Hamil- 
tonian H{u) (j27j) via the two-qubit gate is obtained to be time-dependent of 
the parameter u after some algebra using the formula (j24|) . but the Hamiltonian 
via the the formalism (|22p of this gate is found to be time-independent Heisenberg 
interaction. 

Consider the other solution of the Yang-Baxter equation ()19p , 



R{x) 






V -(1-^) 





1 + a; 
\-x 






-{\-x) 
1 + x 




1 -X \ 



1+x / 



(25) 



where the unitarity condition (j20p requires x real with the normalization factor 
/9 = 2(l + x^). In terms of the variable Q defined by 



x 



cos ( 



sm ( 



X 



2' 



(26) 



the two-qubit gate R(x) ()25p has the formalism of 0, and the Hamiltonian 
(l27ll is calculated to be 



id, 1 , 1 i 1 



cr,. 



(27) 



with the conventional form of the Pauli matrices. The time evolutional operator 

as well as the two-qubit gate ([25]) is given by 



has the form V{Q^ = e 

m 

At 6* = 0, the two-qubit gate R{0) given by 



cos(- 



+ 2isin(--0)/7 = e*(t-2^)^. 
4 



(28) 



^(0) 



/ 1/V2 1/V2 \ 

-1/V2 

1/V2 1/V2 

V -1/V2 1/V2 J 



(29) 



is a unitary basis transformation matrix from the product base to the Bell states. 

Besides the R{x) matrix (j25p . in [7|, the following type of nontrivial unitary 
solutions R{x) of the Yang-Baxter equation (I19p . 



R{x) 



( uxix) uj-jix) \ 

W5(x) CJ3(x) 

LO/^ix) L>Jio{x) 

V a;8(x) W2(x) / 



(30) 
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have been recognized as two-qubit gates modulo a phase factor as well as the 
associated time-dependent Hamiltonians H{x) (j24p also respectively calculated. 
The two-qubit gate R{x) ([30]) is determined by two submatrices Ui and U2, 

U,{x) = ( ''f] ] , U,ix) = ( ''f ] ^f] ) , (31) 

\y UJs{x) U}2{X) J V ^^4(2;) U)q{x) J 

which are in the unitary group U{2). 

According to Terhal and DiVincenzo's understanding [12] of Valient's work on 
matchgates, when the above unitary matrices Ui and U2 are elements of the spe- 
cial unitary group SU{2), quantum computations with the R{x) gate ([30|) acting 
on nearest-neighbor qubits can be efficiently simulated on a classical computer, 
whereas the R{x) gate (|30p combined with the Swap gate acting on farther- 
neighbor qubits may be capable of performing universal quantum computation. 
Obviously, the two-qubit gate ()25p has its two submatrices Ui and C/2 in the 
SU{2) group. 

In terms of Ui and U2 in SU{2), the other two-qubit gate with six nonvanishing 
entries given by 

\ 

■ 

sinh(7 + iO) / 

with the normalization factor p = sinh^ 7 + sin^ 9, can be found in [7] as a unitary 
solution of the Yang-Baxter equation ([9]) with the spectral parameter 9 and real 
parameter 7, but which gives rise to a time-dependent Hamiltonian (j27p . 

However, the two-qubit gate (|22p associated with the delta-function interac- 
tion ([1]) or Heisenberg interaction consists of the following submatrices Ui and 
U2 given by 

U, = -(l U2 = -e--( (33) 

\ 1 y \ I simp cos if J 

where Ui is in SU{2) but U2 is in U{2) because of the global phase. With the 
encoded logical qubit, it has been numerically verified [IT] that the two-qubit 
gate (j22]) itself can achieve universal quantum computation by acting on nearest- 
neighbor qubits. 

4 Concluding remark 

Feynman is well known in the community of quantum information science due to 
his pioneering work on both universal quantum simulation in 1982 and quantum 
circuit model in 1986. Shortly before he passed away in early 1988, very unex- 
pectedly, Feynman wrote: "I got really fascinated by these (l+l)-dimensional 



/ sinh(7 — i9)) 











e*^ sinh 7 
~i sin 9 






—i sin 

-ie 



sinh 7 
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models that are solved by the Bethe ansatz", see Batchelor's feature article [13] 
on the research history of the Bethe ansatz. 

A definition of integrable quantum computation in both physics and mathe- 
matics has been proposed in this article, in order to remove potential conceptual 
confusions between the papers [3] and [7j as well as declare integrable quantum 
computation as an independent research subject in quantum information and 
computation. Note that integrable quantum computation has been also argued 
from the viewpoint of the Hamiltonian formalism of quantum error correction 
codes [Ej. 

The asymptotic condition R{x = 0) of the solution of the Yang-Baxter equa- 
tion (|19p satisfies the braid group relation and hence R{x = 0) can be viewed 
as a unitary braiding gate, which suggests similarities and comparisons between 
integrable quantum computation and quantum computing via unitary braid rep- 
resentations, see [71 [U] for the detail and [3j for comments. 
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